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ABSTRACT 

A theoretical model is developed for making computerized fore- 
casts of mixed layer depth. An empirical equation relating significant 
wave height to layer depth is used to account for mechanical mixing. 
Heat exchange parameters are introduced into the model to produce 
convective mixing. The u and v components of surface current aTe 
used to investigate the effect of convergence and divergence. The 
layer depths forecast by the model are compared to actual layer depths 
and to the Fleet Numerical Weather Central layer depth analysis. 
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I. INTRODUCTION 



With the increased emphasis on antisubmarine warfare and the 
development of various sonar range prediction models, it has become 
increasingly important that an accurate prediction of the ocean thermal 
structure be made available. Changes in ocean thermal structure, 

-V* 

which occur mainly in the surface layer and the thermocline, are of 
major importance in determining variations in the vertical sound 
velocity profile. The sound velocity field determines the refraction 
pattern and the existence of surface ducts, sound channels and other 
sound transmission paths. Therefore, better data on the thermal 
structure will increase the reliability of predicting sonar performance. 

The most important thermal structure parameter is the vertical 
extent of the surface duct, or the mixed layer depth. A surface duct 
will exist if: (1) the temperature increases with depth, or (2) an 

isothermal layer extends downward from the sea surface. In case 
(1), sound velocity increases with depth as temperature and pressure 
increase; in case (2), there is no temperature gradient and, in the 
absence of a salinity gradient, pressure causes sound velocity to 
increase with depth. The greater the depth of the mixed layer, the 
greater is the difference between the sound velocity at the surface and 
at depth. The number of sound rays which are trapped in the layer is 
therefore greater. Also, the deeper the layer, the fewer the surface 
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reflections which are required to reach a given range, with a resultant 
decrease in propagation losses. With increasing layer depth, the 
detection range increases for a target at its best depth to avoid 
detection. Thus it is evident that the mixed layer depth, to a great 
extent, determines the sonar ranges which will prevail. These 
variations in ray propagation with layer depth are illustrated in 
figure 1. 

The mixed layer depth forecasting scheme presently in use at 
Fleet Numerical Weather Central (FNWC) is based primarily on 
climatology and analysis of previously reported conditions. Each day, 
the previous forecast field is adjusted towards climatology and then is 
modified by turbulent mixing due to wave action. This method of 
predicting mixed layer depth is rather elementary: the use of this 

model gives very little improvement over a prediction based on 
climatology or persistence (forecast conditions will be the same as 
present conditions). An exception occurs when wind mixing causes a 
significant change from climatology; under these conditions a fair 
correlation exists between the forecast and actual conditions. A more 
complete explanation of the method presently in use at FNW C is con- 
tained in Appendix A. 

The purpose of this research is to produce a computerized mixed 
layer depth model for FNWC which will make a prediction of layer 
depth based on environmental conditions to replace the present one 
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Figure 1. Variation of sound ray propagation with layer depth. 
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which is strongly dependent on climatology. The prediction method 
will be considered improved if the mixed layer depth produced proves 
to have a smaller percent error when compared to the verification 
layer depth than the layer depth produced by the FNW C analysis. A 
further check as to the accuracy of the prediction can be made by com- 
paring the standard deviation of the forecast layer depth from the 
verification layer depth to the standard deviation of the FNWC method: 
25 feet in the summer and 80 feet in the winter. 
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II. PREVIOUS MIXED LAYER DEPTH STUDIES 



Numerous studies relating to various aspects of the mixed layer 
depth have been conducted at the U. S. Naval Postgraduate School. 

Clark [4] directed her attention to the upper layers of the atmos- 
phere and such conditions that might possibly be associated with 
variations in the temperature structure of the ocean. She presented 
two hypotheses: (1) the upper mixed layer of the ocean and the layer 

of air in immediate contact with the ocean surface should be treated 
as a unit which is affected by meteorological parameters above this 
unit, (2) there exists in the upper atmosphere a ’mirror image’ level 
where fluctuations in wind speed closely depict oscillations in the depth 
of the mixed layer. Clark reached the conclusion that a 'mirror image' 
level does exist in the atmosphere at an altitude of 6000 meters. 

Geary [10] developed a model of mixing by wind generated wave 
motion and compared it to Laevastu’s [ 1 8] and Neumann's [23] methods. 
He concluded that Laevastu’s method gave the closest approximation to 
an upper limit on the observed mixing. He also stated the need for the 
development of a convective mixing model. 

Edgren and MacPherson [8] used multiple regression analysis on 
various parameters (primarily meteorological) to examine the physical 
processes causing increases and fluctuations in the mixed layer depth. 
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It was determined that convection is the process which causes the 
seasonal increase of the layer depth during the cooling season and that 
short term fluctuations of the mixed layer are due primarily to internal 
waves . 

Luskin [2 1] d eveloped a simple model for mixed layer depth change 
by thermohaline convection and outlined a method for adapting the con- 
vection model for use as a forecast tool. He pointed out the need to 
develop a heat budget model for forecasting convective mixing. 

Lambright [20] investigated the apparent random oscillation of the 
mixed layer depth and by the nature of its energy spectrum and its 
correlation with tidal activity, determined the source of some com- 
ponents of this oscillation. 

Davis [5] investigated Bulgakov's model relating the depth of con- 
vective mixing to density and salinity changes in the upper layer of the 
ocean. He suggests the use of this model in taking into account the 
convective mixing due to increased density caused by mixing two or 
more layers of the ocean during the heating season. 

Hancock- [12] verified the Fleet Numerical Weather Central 
analyses and forecasts of mixed layer depth. The model verified in 
this study included convergence/divergence computations. The con- 
vective mixing influence was calculated from a change in sea surface 
temperature. 

Ciboci [3] verified Mazeika's [22] method for Ocean Weather 
Station "P". It was found to be successful during the heating season 
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followed by a rapid decline as the cooling season began. He also 
determined that the forecasting curves developed by James [15] are 
easier to use than Mazeika's and provide a more accurate mixed layer 
depth forecast. 

Kelley [ 1 6 ] studied a model based on Kitaigorodsky 1 s application 
of similarity theory and modified by McDonnell to forecast mixed 
layer depth* 

In this present work, an attempt is made to use various tech- 
niques previously developed to produce a forecasting method which 
takes into account heat budget calculations in computing the effect of 
convective mixing. The need for a model of this type has been stated 
many times and this study is a further step in the continuing quest to 
develop an accurate forecasting model for mixed layer depth. 
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III. FACTORS AFFECTING MLD 



The ocean may be considered to consist of three distinct layers: 

(1) the mixed upper layer, influenced by atmospheric conditions; 

(2) a zone of rapidly decreasing temperature called the thermocline 
also influenced by atmospheric conditions, and (3) a region of slower 
changes extending from the thermocline to the ocean bottom. 

It is a feature of the mixed upper layer, specifically the depth at 
which this layer ends and the thermocline begins, that is of primary 
interest in this work. This depth is termed the mixed layer depth (or 
simply layer depth) and is defined as the depth of maximum tem- 
perature or the extent of the isothermal layer. The term isothermal 
is defined to include all positive gradients, and negative gradients up 
to and including 0. 3 F/100 feet [13]. This definition is equivalent to 
the definition of sonic layer depth, the depth at which maximum sound 
velocity occurs in the layer above the thermocline [2]. Thus the two 
terms may be used interchangeably. 

A number of methods are available for predicting the thermal 
structure of the oceans. These involve calculations of net heat ex- 
change across the air- sea boundary, convective and wind mixing, 
continuity, and advective processes. Each procedure may be utilized 
separately, but the most accurate forecasts are generally made by 
considering the collective effects of all processes. 
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A. NET HEAT EXCHANGE 



The sea- surface temperature and the thermal structure imme- 
diately below the surface and its variations are greatly influenced by 
the heat exchange between the sea and the atmosphere. It is this energy 
exchange at the air-sea interface which forms the basis for the heat 
budget method of forecasting near- surface thermal conditions. The 
net heat gain or loss at the ocean’s surface for the period of the fore- 
cast is computed from known or forecast values of sun's altitude, 
cloud cover, humidity, sea- surface temperature, wind speed, relative 
vapor pressure, and albedo. Of these parameters, only the sun's 
altitude is a truly-known quantity, while the others are forecast with 
varying degrees of accuracy. Heat budget computations are also of 
use in solving other problems, such as ice prediction, prediction of 
transient thermoclines and possibly the improvement of medium- range 
weather forecasts over the sea. 

Once the total heat exchange has been calculated, the gain or loss 
of energy must be distributed throughout the water column. Qual- 
itatively, heat gains often are associated with decreased layer depth 
while losses accompany layer depth increases. 

B. MIXING 

The primary types of mixing considered in the formation of the 
mixed layer are: (1) instability (or convective) mixing produced by 

the sinking of dense water; (2) mixing as the result of convergence or 
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divergence causing the sinking or rising of water; and (3) mechanical 
(or forced) mixing, a turbulent transfer of momentum from one level 
to another by the combined action of wind waves and associated wind 
currents. All three types of mixing often occur simultaneously. 

After calculating the magnitude of the heat gain or loss, convective 
mixing and mechanical mixing must be considered to determine heat 
distribution. For a net heat gain, implying a negative temperature 
gradient in the surface layer, wind mixing becomes the dominant 
factor in computing the mixed layer depth. 

Although wind is the direct force involved in mechanical mixing, 
mixing in the water column is accomplished through the action of 
secondary phenomena- wind waves and pure wind currents. Swell 
can also contribute to mechanical mixing when combined with the 
effects of wind waves and currents. The mixing by waves is due 
primarily to the presence of waves of different period and height and 
to the turbulence caused by breaking. In the majority of cases, the 
seas are not fully developed, and wind speed alone cannot be used 
directly in the determination of mechanical mixing by the sea. In 
determining the thickness of a mixed layer created by wave action, 
the surface wave characteristics are first determined and then a 
depth is computed where particle motion and the resulting mixing are 
negligible. 

The mixed layer depth due to mechanical mixing will depend on 
the wave model used. Since the singular wave forecasting model 
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(one which relates wave height to wind speed through a simple exponen- 
tial relationship) is presently in use, only the simplest relationship 
between wave height and mixed layer depth is justified [19]. 

C. CONVECTIVE MIXING 

Although mechanical mixing is an important parameter in deter- 
mining the mixed layer depth, it is not the only one. In fact, it may 
at times be entirely masked by other factors, such as convective 
mixing and convergence/divergence. 

When the net heat exchange is from ocean to atmosphere (negative 
heat exchange), convective mixing is the dominant factor in determining 
the mixed layer depth. The upper layers of the ocean become unstable 
due to cooling and sink until they reach a stable level. This convection 
will result in isothermal mixing of the water to the level of stability. 

After determining the effects of convective mixing, the wind- 
mixed layer depth is computed. If wind mixing is deeper than convec- 
tion, the column is mixed isothermally to the wind-mixed depth. If 

\ 

the convective layer is deeper, no further modification is required. 

D. CONVERGENCE/DIVERGENCE 

Pure wind currents may exert another influence on the layer in 
addition to the part they play in forced mixing. In an area in which 
there is a wind-created divergence of surface water (waters horizon- 
tally advected apart) the layer depth decreases since cooler water 
from intermediate depths must well upward to replace the divergent 
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surface waters. Similarly, a wind- produced convergence (waters 
flowing together) will cause an increase in the layer depth due to the 
sinking of the accumulated surface water. 

E. ADVECTION 

Advective changes to the thermal structure are the most difficult 
to define with any degree of accuracy. The degree to which advection 
will affect the thermal structure is dependent upon the length of the 
forecast period, the local circulation, the horizontal temperature 
distribution, and the wind conditions. In areas where advection is 
determined to be of consequence, primary consideration is given to 
drift currents since their effects tend to be dominant. Permanent 
currents are assumed to be geostrophic (determined essentially by 
the temperature distribution if there is negligible salinity contribution) 
and thus produce no heat advection. 

F. INTERNAL WAVES 

A primary factor which limits the predictability of the instanta- 
neous mixed layer depth and verification of the forecast is the presence 
of internal waves on the interface between the mixed layer and the 
thermocline. Characteristically, these waves have large vertical 
displacements at the bottom of the mixed layer. While similar to sur- 
face waves, internal waves lack the turbulent interchange of momentum 
(a result of wave breaking) possessed by surface waves. Thus physical 
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mixing due to internal waves is not great unless these become break- 
ing waves. Therefore, internal waves generally cause a periodic rise 
and fall of the mixed layer depth but do not seriously affect the average 
position of the bottom of the mixed layer. For this reason, as well as 
the lack of sufficient knowledge concerning the spectral distribution of 
internal waves, they were not considered in this study. 
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IV. PROCEDURE 



A. DATA 

In the development of this mixed layer depth forecast procedure, 
it was decided that available bathythermograph (BT) reports would be 
used to obtain the actual mixed layer depth (MLD) for use as a verifi- 
cation of the predicted value. For this reason, the investigation is 
concentrated in the vicinity of permanent ocean station vessels on the 
supposition that these locations would have the highest density of BT 
data available. 

It was also desired to check the forecast model during different 
seasons; the study is confined to data from the months of January and 
September since these are quasi-stationary periods when mean physical 
parameters are relatively invariant. 

First, a two degree latitude-longitude square around each of 
eight ocean station vessels was checked for the availability of BT data 
for ten-day periods in January 1972 and September 1971. Due to data 
limitations, the areas around Ocean Station Vessels (OSV) M E M , "N n , 
and n P n were selected as the locations for use in developing the pro- 
cedure. The periods used were 10 to 20 September 1971 and 10 to 20 
January 1972. Since these positions were not on a grid point for the 
standard FNW C analysis grid network, the data fields were inter- 
polated to the OSV locations using linear interpolation between the 
four surrounding grid points. 
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For the periods selected, the FNW C fields of analyzed sea surface 
temperature, forecasted total heat exchange for the next 24 hours, 
u and v components of surface current, analyzed sea temperature at 
600 feet, and forecasted wave and swell height were obtained. These 
fields were the input data to the forecast procedure. A brief descrip- 
tion of each field is contained in Appendices B through F. 

The order followed in producing the forecast can be seen in figure 
2. The computer programs developed to produce the MLD forecast are 
included as Appendix H. 

The first step in the process is to input an initial layer depth. 

This depth will be extracted from a BT taken 24 hours prior to the 
desired forecast time. The next step in the initialization phase is to 
input the data to be used in the forecast computation. 

B. FORCED MIXING (HEAT EXCHANGE POSITIVE) 

The actual forecast procedure starts with the determination of 
whether the net heat exchange for the next 24 hours is positive or neg- 
ative. If the net heat exchange is positive, the layer depth due to forced 
wind mixing is computed. Three methods of determining the depth of 
turbulent mixing are tested to determine which will give the best results. 

The first method is an empirical formula based entirely on the 

significant wave height (H-^^h suggested by Laevastu [ 19 ].* 

MLD = 12.5 H (1) 

The significant wave height used in this formula is actually the com- 
bined wave height (CH) due to both wave height (WH) and swell height 
(SH), calculated as follows: 
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Figure 2. Mixed layer depth model flow diagram. 



( 2 ) 



CH =WH 2 + SH 2 

The second method of forced mixing investigated is the one 
presently in use at FNWC. This also is based primarily on the com- 
bined wave height (CH), with a slight correction for the magnitude of 
the temperature gradient between the surface and 600 feet [9]: 



MLD = CH [ 11 - 0. 1 (SS T - T 6oo ) J 


SST - t 600 


MLD = 11- CH 


SST < T^oo 



The final scheme tested is that suggested by James [15]. The 

first step in this method is to calculate the absorption of radiation by 

layers. Absorption coefficients by layers for various water types 

(clear oceanic, average oceanic, average coastal) are given by James. 

For the purposes of this study, average oceanic turbidity is assumed. 

This is a valid assumption since this model is not intended for use 

in coastal waters. For oceanic water, with a net heat exchange at 

-2 

the surface of 500 cal-cm , the resulting temperature difference 

between clear oceanic and average oceanic water in the 0-20 foot 
o 

layer is 0. 08 F. This difference is not of sufficient magnitude to 

justify the additional computer time required to make a determination 

of water type. In addition, absorption will be computed for only the 

0-20 foot layer. This is done to decrease the number of computations 

needed and is considered valid, since an input at the surface of 
r - 2 

500 cal-cm is required to produce a change in temperature of 
o 

0.1 F in the 20-40 foot layer while this same energy input at the 
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o 

surface would cause a change of 1. 3 F in the 0-20 foot layer. It is 
the gradient in this shallow layer which must be overcome by mechanical 
mixing. The average change in temperature in the layer due to the 
surface heating is calculated, and is then added to the original tem- 
perature at the ten-foot level. A linear gradient is assumed through 
this point between the surface and the original trace at 20 feet. 

The next step is to compute a mixing parameter (M ) which is a 
function of significant wave height and the period of maximum 

energy (T max ) [ 15]: 

•^w = ^max^l/3 * W 

The T max in this case is the period of maximum energy which would 

be present for a fully developed sea having the given significant wave 
height. The combined wave height (equation 2) is again used in this 
case as the significant wave height. Once this mixing parameter has 
been computed, it is used with the temperature gradient of the 0-20 
foot layer to compute the layer depth [ 15 ] : 



LD = 



Kl 

a T 



( 1 - e 



-K? aTM 

4 W\ 



where and are constants. 



C. CONVECTIVE MIXING 



( 5 ) 



For the situation of a net heat loss to the atmosphere, layer 
depth change due to convective mixing must be calculated. This is 
accomplished using a method set forth by James [15]. The first step 
is to calculate the convective mixing parameter, M^: 
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M = 



2 Q! 



C C p?w - T 



( 6 ) 



-2 



where Q is the net heat loss in gram cal. -cm , C is the specific 
1 P 

o 

heat of sea water and is assumed to be constant at 0. 935 cal/gm/ C, 

( 3 ^ is the density of sea water and is assumed to be constant at 

3 

1.024 gm/cm , and aT is the average temperature gradient from 

the bottom of the initial layer to 600 feet in F/100 feet. M is then 

c 

combined with the initial layer depth (h^) to give the mixed layer 
depth due to convection as follows: 



LD = AAV 2 + (M C )2 * (7) 

The next step is to calculate the forced mixed layer depth using 
one of the methods described in the previous section. The mechanical 
and convective mixed layer depths are then compared and the deeper 
one is taken to be the actual layer depth. In various references, it 
has been suggested that when convective mixing is present, mechanical 
mixing can be ignored. This suggestion is also tested. 



D. CONVERGENCE/DIVERGENCE 

The effect of convergence or divergence on the layer depth is 
then determined. This is accomplished through the use of 
the u and v components of surface current, which are combined to 
give a change in layer depth due to convergence or divergence as 
follows [12]: 3 4 

I 2 . 
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LD = 0. 864 (U 1 + U 3 -U 2 -U 4 +V 3 +V 4 -V 1 -V 2 ) LD/4L 
where 0. 864 is a dimensional constant to change the variation in the 
layer depth from the metric system to feet, and L is the distance 
between grid points. The decision to work in feet rather than meters 
was made in order to decrease the number of calculations (wave 
heights, swell heights, layer depths and temperature gradients are 
all in feet). This value for the change in layer depth is then added to 
the layer depth previously calculated and the result is the predicted 
mixed layer depth. 

E. ADVECTION 

The advection of heat and the change in the thermal structure due 
to this advection is then computed. 

For the areas under investigation, it was noted that due to the 

o 

slight temperature gradient between grid points (maximum of 0. 01 F 

per nautical mile) and the small currents (maximum transport on the 

order of 50 nautical miles per day) the temperature change due to 

o 

advection would be no more than 0. 2 F/day at any time. 

Due to the difficulty in describing advection accurately and the var- 
iable contribution it makes in the areas under investigation, the 
decision was made to disregard advection in this study. 

F. FITTING THE FORECAST TO DATA 

At this point in an operational procedure, the latest BT reports 
are read into the program and the forecast field is adjusted 
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accordingly with respect to the reported layer depth and the age of 
the BT report. In this study, this step was omitted since the actual 
BT reports were the basis for determining the accuracy of the 
forecast. 
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V. RESULTS 



A. COMPARISON OF MECHANICAL MIXING FORMULAS 

In order to determine the mechanical mixing method to be used in 
the procedure for the prediction of mixed layer depth, calculations of 
forced mixing were made using each of the three proposed relation- 
ships. In relat 1 to the complete forecasting model (figure 2), the 
calculations in this section refer only to the forced mixing step. The 
heat exchange is given as positive. The data input are combined wave 
height and the temperature gradient between the surface and 600 feet. 
All remaining steps in the model are ignored. The computations were 
made for five specified periods when mechanical mixing was the 
determining factor in the forecast layer depth. This is the case when 
the heat exchange is positive. The combined wave height (CH) was 
also a factor in selecting the period to be investigated; computations 
were desired for a range of wave heights. The layer depth from each 
relationship and the observed layer depth are compared in Table I. 

The percent error from the actual conditions (in parentheses) and the 
combined wave height for the forecast period are also shown. 

For reference, the equations of the three relationships are as 
follows : 

1. Laevastu 

LD = 12. 5 CH (9) 
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2 . 



FNWC 



LD = CH [11- 0. 1(SST- T 6qo )] 



SST " T 600 < 3 > 



LD = 11- CH 



SST*= T^qq 



3. James 




(5) 



AT 



A comparison of the Laevastu formula to the FNWC one shows that 
the FNWC will always produce a shallower forecast. The FNWC con- 
stant (11) is 12% smaller than the Laevastu constant (12.5) thus pro- 
ducing a forecast layer depth due to mechanical mixing which, in the 
absence of a temperature gradient between the surface and 600 feet, 
will be 12% less than the forecast produced by the Laevastu formula; 
as the temperature gradient between the surface and 600 feet increases 
in magnitude, the differences between the layer depths forecast by the 
two methods increases. 

In comparing the Laevastu results to those of James it was noted 
that the percent difference between the two layer depths varied with the 
wave height. With wave heights of the order of five feet, the difference 
between the two models was less than 10% (absolute differences in the 
range 1-5 feet). As the wave height increased to 10-20 feet, the James 
model produced layer depths approximately 18% greater. With a large 
wave height (26. 6 feet) the Laevastu model produced a layer over 50% 
deeper than the James model, and 18% shallower than the observed 
layer depth. 
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TABLE I. 



Comparison of layer depths in feet produced by various 
mechanical mixing formulas. Numbers in parentheses 
are percent error from observed conditions. Combined 
wave height (CH) in feet and net heat exchange (Q ) in 
cal/cm are indicated for each case. is positive in 

all cases . 





CH 


Q 

n 


LAEVASTU 


FNWC 


JAMES 


OBSERVED 


1 . 


4. 3 


281 


54 (51. 9) 


41 (62. 7) 


50 (54.5) 


110 


2. 


4. 7 


67 


59 ( 1.7) 


45 (25. 0) 


60 ( 0 ) 


60 


3. 


10.4 


25 


131 (69.5) 


112 (74.0) 


160 (62. 8) 


430 


4. 


18.5 


257 


231 (43. 5) 


200 (51. 1) 


270 (34.0) 


409 


5. 


26. 6 


308 


332 (18. 2) 


287 (29. 3) 


160 (60. 6) 


406 
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One aspect of the three mechanical mixing relationships which 
should be noted is that as the net heat gain at the ocean’s surface 
becomes larger, the FNW C and James models produce shallower layer 
depths due to the influence of the increased thermal gradient in the 
upper layer. The Laevastu formula contains no consideration of heat 
gain. 

Only in case 2 (Table I) was an accurate prediction of layer depth 
under the influence of mechanical mixing produced (by the Laevastu 
and James relations). Cases 1 and 2 occurred in September during a 
period of decreasing surface heating and low wave height when the 
layer depths were relatively shallow (60 to 160 feet). In these cases, 
the wind mixing was not sufficient to keep the surface layer mixed 
below 60 feet and after two days of surface heating, a shallow transient 
layer was formed (the layer at 60 feet was destroyed the following day 
by convective mixing). 

The Laevastu relation exhibited fair results (18% error) in case 5. 
Cases 3 through 5 occurred in January during a period of increasing 
surface heating and large combined wave height when the layer depth 
was fairly deep (400 feet). For these cases, although the surface 
heating was large, the mechanical mixing was of such magnitude that 
a shallow surface layer could not form. 

It appears that there is some relationship between mechanical 
mixing, heat exchange, and the initial layer depth which determines 
whether a layer will be formed which conforms to the mechanical 
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mixing term or if the mechanical mixing will act to maintain the 
original deep layer by preventing the formation of a shallow layer. 
When the conditions are such that the mechanical mixing acts to 
eliminate the shallow layer, none of the mechanical mixing methods 
investigated here produce accurate forecasts for the deeper layer. 

After a consideration of the preceding factors, it was decided 
to use the Laevastu formula for predicting mechanical mixed layer 
depth in the model tested in this study. The Laevastu formula was 
chosen over the FNW C on the basis of the smaller error between the 
forecast and observed layer depth. The James formula was eliminated 
primarily due to the difficulty in computerizing it (computations in 
Table I were produced manually). This was not the only considera- 
tion however, as the accuracy of the two methods (Laevastu and 
James) was also compared. At wave heights up to 20 feet, there was 
little difference in the two while in the fifth case (wave height of 
26. 6 feet) the Laevastu result was much better (18% vs. 60%). 

The difficulty in computerizing the James model arises from the 
fact that the constants (K^ and K^) in equation 5 vary with the tem- 
perature gradient. One possible solution for this is to use a nomo- 
gram provided by James [ 1 5 ] to produce a grid of layer depths with 
the wind mixing parameter as one ordinate and the temperature 
gradient as the other. This could be done but the computer time 
required for the grid search in each case would be prohibitive. 
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B. DAILY FORECASTS 



Using the various procedures outlined in Section IV, mixed layer 
depth forecasts were produced for selected 24-hour periods. In 
making the forecasts in this and the following section, each applicable 
step of the model (figure 2) was considered (the path taken depending 
on whether the heat exchange was positive or negative). For the 
forecasts in this section, the initial layer depth was determined from 
the BT report taken 24 hours prior to the forecast time. The criterion 
for the selection of a forecast period was the availability of BT data 
for use as the initial layer depth and as a verification of the accuracy 
of the forecast. The forecasts are shown in the computer output sec- 
tion (Appendix I). 

The purpose of this phase was to check the accuracy of the fore- 
cast over a 24-hour period using an accurate layer depth as an input. 
The observed layer depth, forecast layer depth and the FNW C analysis 
of the layer depth at the forecast time for each forecast are shown in 
Tables II thru V. The percent error between the forecast and observed 
layer depth and the FNWC analysis and the observed layer depth were 
calculated for each case and are displayed as the value in paren- 
theses in each table. The last line of each table shows the average 
standard deviation of the forecast layer depth from the observed 
layer depth and the FNWC analyzed layer depth from the observed 
layer depth. 
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1 . January, Station n N n 



The results for this forecast period are contained in Table II. 

These forecasts were made for *a period of substantial negative heat 

2 2 

exchange (-400 cal/cm to -200 cal/cm ) with combined wave heights 
of 6-11 feet. The procedure exhibited improved accuracy in two of the 
four cases. For the periods when the FNW C analysis was more ac- 
curate, the improvement was only about 2% in each case. The max- 
imum error in the procedure was 7. 9% while the analysis had two values 
with an error greater than 8%. With respect to predicting the average 
conditions for the period, the procedure exhibited an error of 3.6% 
with an average standard deviation of 29 feet as compared to -7. 5% 
and 53 feet for the FNWC analysis. 

The FNWC analysis was less than the observed conditions 
in every case while the procedure forecast was higher in all but one 
case. 

2, January, Station ,, P" 

Results are shown in Table III. The first two forecasts 

2 2 

were for periods of small heat exchange (-90 cal/cm and 25 cal/cm 

respectively) while the remaining four were for a period of sub- 

2 

stantial negative heat exchange (-250 cal/cm ). The combined wave 
height for the period was 10-13 feet. For this period, the FNWC 
analysis proved to be more accurate in four of the six cases and close 
to the model (less than 2% difference in percent error) in the remaining 
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TABLE II. 



Comparison of observed MLD, forecast MLD, and 
FNW C MLD analysis (in feet) for Station M N M , January 
(numbers in parentheses are percent error). 



TIME 


DATE 


OBSERVED 


FORECAST 


FNW C 


00 


12 


530 


552( 4.2) 


483 ( — 8. 


12 


12 


5 70 


5 6 1 ( - 1.9) 


499(- 12. 


00 


13 


510 


53 1( 4.1) 


502(_ 1. 


12 


13 


530 


572( 7.9) 


497(- 6. 


AVERAGE 




535 


554( 3.6) 


495(- 7. 


STD. DEV 






29 


53 
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TABLE III. Comparison of observed MLD, forecast MLD, and 

FNW C MLD analysis (in feet) for Station M P n , January 
(numbers in parentheses are percent error). 



TIME DATE 


OBSERVED 


FORECAST 


FNW C 


00 


13 


360 


400( 11. 1) 


3 9 1 ( 


8. 6) 


*12 


13 


430 


1 3 0 ( 69. 8) 


404( - 


6. 0) 


00 


39 


390 


440( 12.8) 


389(- 


0. 2) 


12 


19 


410 


3 86 ( - 5.9) 


404(- 


1.5) 


00 


20 


380 


403 ( 6. 1) 


405( 


6. 6) 


12 


20 


460 


420(- 8.7) 


4 12( _ 


10. 4) 


AVERAGE 




405 


363( - 10. 4) 


400(- 


1. 2) 


STD. DEV. 






139 


30 





* indicates positive heat exchange 
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two cases. The maximum error for the FNWC analysis was 10.4% 
while that for the procedure was -69. 8%. The procedure failed to 
produce a forecast which was within 5% of the observed layer depth 
for any case. 

In forecasting the average conditions, the procedure had an 
error of -10.4% with an average standard deviation of 139 feet while 
the FNWC analysis had an error of only -1.2% and an average standard 
deviation of 30 feet. 

3. September, Station n N ,! 

Results are shown in Table IV. The net heat exchange for 

the first two forecasts is positive; the heat exchange for the remainder 

2 2 

of the period is negative, increasing from -83 cal/cm to -200 cal/cm 
during the period. The combined wave height throughout the forecast 
period is 2-6 feet. In comparing the two schemes, the results of 00Z 
on the 19th are discounted since, although the FNWC analysis had a 
smaller error, neither was within 500% of the observed layer depth 
(this value was included in the calculation of the average layer depth 
and the average standard deviation). Of the remaining 10 forecasts, 
the FNWC analysis proved to be more accurate in six cases; however, 
the best result obtained was an error of 7. 9%. 

The procedure showed an error of 15% in forecasting average 
layer depth with an average standard deviation of 85 feet. The FNWC 
analysis had a 3. 8% error in predicting average conditions with an 
average standard deviation of 52 feet. 
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TABLE IV. 



Comparison of observed MLD, forecast MLD, and 
FNW C MLD analysis (in feet) for Station "N", September 
(numbers in parentheses are percent error). 



TIME DATE 


OBSERVED 


FORECAST 


FNW C 


*00 


13 


110 


53( 


- 51.8) 


132( 20. 0) 


* 12 


13 


60 


59( 


- 1. 7) 


143(138. 0) 


00 


14 


180 


1 10( 


- 38.9) 


141(-2 1. 7) 


12 


14 


160 


60( 


- 62.5) 


141( 9.2) 


00 


16 


150 


1 5 1 ( 


0. 7) 


1 3 8 ( - 8.0) 


12 


16 


180 


132(. 


- 26. 7) 


138(-32. 3) 


00 


19 


20 


1 9 1 ( 


855. 0) • 


135(575. 0) 


12 


19 


140 


160( 


14. 3) 


129(- 7.9) 


00 


20 


170 


28( 


83.5) 


137(-25. 3) 


12 


20 


170 


141( 


17. 1) 


137(- 19. 4) 


AVERAGE 




133 


1 1 3 ( 


- 15.0) 


138( 3.8) 


STD. DEV. 








85 


52 



* indicates positive heat exchange 
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4. 



January, Station n E M 

Results are presented in Table V. The net heat exchange for 

2 

the first two forecasts is approximately -400 cal/cm ; for the remain- 

2 

ing forecasts the heat exchange is approximately -70 cal/cm with the 

exception of 00Z on the 16th when there was a positive heat exchange 

2 

of 9 cal/cm . The combined wave height during the period is 3-7 feet. 
For the first forecast layer depth, both methods show a large error; 
in the remaining five cases, the procedure presented in this study is 
more accurate four times. 

In predicting average conditions, all forecasts were included 
(the forecast for 00Z on the 12th was not considered in computing the 
average standard deviation). The procedure exhibited an error of 
36. 7% with an average standard deviation of 80 feet while the FNW C 
analysis had an error of 6% and an average standard deviation of 
240 feet. 

5. Summary 

Large deviations of the results from observed conditions 
occur primarily when the net heat exchange at the ocean surface 
becomes positive for a short period of time. In this case, the layer 
depth forecast by the model is determined entirely by mechanical 
mixing. When this happens, the result is the prediction of the forma- 
tion of a shallow layer (as the 12Z forecast on the 13th in Table III 
and the 00Z forecast on the 20th in Table IV show) but this may not 
occur. To correct this problem, some method must be devised to 
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TABLE V. 



Comparison of observed MLD, forecast MLD, and 
FNWC MLD analysis (in feet) for Station M E n , January 
(numbers in parentheses are percent error). 



TIME 


DATE 


OBSERVED 


FORECAST 


FNWC 


00 


12 


110 


910(727. 0) 


422(284. 0) 


12 


12 


580 


6 6 1 ( 14.0) 


434( -25 . 2) 


00 


15 


620 


770( 24.2) 


5 3 1 (- 14. 4) 


12 


15 


650 


679( 4.5) 


544(- 16. 3) 


*00 


16 


100 


47 ( 53.0) 


493(393. 0) 


12 


16 


670 


666(- 0. 5) 


464(-30. 7) 


AVERAGE 




455 


622( 36. 7) 


4 8 1 ( 5.7) 


STD. DEV. 






80 


240 



* indicates positive heat exchange 
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determine whether this predicted shallow layer will actually form or 
whether the deeper layer will persist. This could perhaps be related 
to some value of wind mixing which will cause turbulence to such an 
extent that a shallow layer cannot form. Another method would be 
to set a limit (either absolute or a percentage) on the decrease in 
layer depth which would be permitted to occur. When this limit was 
exceeded, the forecast mixed layer depth would be forecast by persist- 
ence (the forecast layer is the same as the present layer depth). 

Large errors between the FNWC analysis and the observed 
conditions occurred due to the deviation of the observed conditions 
from climatology. The variability of the climatological data is 
presently being evaluated by FNWC and others. 

Although internal waves were disregarded in this study, it is 
of value to note their possible effect on the accuracy of the forecast. 
Assuming an average amplitude of 40 feet for the internal waves 
affecting these forecasts, the forecasts can be considered to be 100% 
correct if they are within +40 feet of the average observed layer depth. 
Applying this assumed internal wave amplitude to the forecasts yields 
an accuracy of 90% in the first two cases (Tables II and III). For the 
other two periods, the improvement is negligible, indicating that 
perhaps other factors are present or that larger internal waves are 
pr es ent. 

Examination of the average standard deviations indicates 
that although the average standard deviation may be large, both the 
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procedure and the FNW C analysis should be able to predict the average 
conditions with some degree of accuracy. 

C. CONTINUOUS FORECASTS 

The purpose of this phase was to see if the procedure could take 
an actual layer depth (from the BT report taken 24 hours prior to the 
forecast time) as a starting point and make successive forecasts for 
a period of time using the forecast mixed layer depth as the starting 
point for each succeeding prediction. It was hoped that after running 
the model for a period of time in this way, the difference between the 
resulting forecast layer depth and the observed condition would be 
minimal. The forecasts are shown in the computer output section 
(Appendix I). ' 

The criterion used in selecting the periods used was the avail- 
ability of observed layer depths throughout the period for use in making 
a comparison of the accuracy of each forecast. The forecast MLD, 
observed MLD, and the percent error between the two are shown in 
Tables VII thru X. 

1. January, Station M N M (00Z forecasts) 

Results are shown in Table VI. The procedure worked well 
(maximum error of 8. 4%) for the first three days. On the fourth day, 
there was a positive heat exchange, thus causing convective mixing to 
cease and the forecast layer depth was determined by mechanical 
mixing. Since mechanical mixing was small on this day, the layer 



45 



TABLE VI. 



Comparison of the forecast MLD to the observed MLD 
(in feet) and the percent error between the two for 
Station "N", January (00Z forecasts). 



DATE 

12 

13 

14 
*15 

16 

17 

*18 

19 

20 

AVERAGE 



FORECAST 

552 

553 

554 
43 
82 

115 

46 

77 

87 

234 



OBSERVED 

530 

510 

520 

540 

440 

540 

520 
500 
590 

521 



% ERROR 

4 . 2 
8 . 4 
6 . 5 
- 92 . 0 
- 81.4 
- 78 . 7 
- 91 . 1 
- 84 . 6 
- 85 . 2 
- 55 . 1 



indicates positive heat exchange 



depth was forecasted to decrease drastically (from 520 feet to 43 feet) 
while actually it increased by 20 feet. Since this layer depth (43 feet) 
was used as the initial layer depth for the next forecast, the large 
error (80%) was propagated throughout the remainder of the period. 
Positive heat exchange again occurred on the 18th which caused the 
layer again to decrease; this stopped the growth which was occurring 
due to convective mixing. 

Calculations were made to determine the effect on the ac- 
curacy of the forecast if for the forecast on the 15th, persistence of 
the layer depth' on the 14th (554 feet) were forecast rather than the 
mechanical mixing layer depth (42 feet). The forecast procedure is 
then followed as before until the 18th when again the layer depth for 
the previous day is forecast to persist. The result of this modified 
forecast is that the procedure predicts the average layer depth with an 
error of 6. 3% and an average standard deviation of 53 feet. These 
values compare favorably with the values computed for the three days 
on which good forecasts were obtained originally (6. 3% and 41 feet). 

2. January, Station f, P n (00Z forecasts) 

Results are shown in Table VII. The results for this period 
were very good although there were three days during the period for 
which no BT observations were available for comparison. 

This period was a time of strong wind mixing (18 to 25 foot 
combined wave height) in addition to being an interval of substantial 
convective mixing. The large mechanical mixing on the 14th and 15th 
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TABLE VII. Comparison of forecast MLD to the observed MLD 
(in feet) and the percent error between the two for 
Station "P", January (00Z forecasts). 



DATE FORECAST 



OBSERVED % ERROR 



13 370 

*14 233 

*15 336 

16 336 

17 340 

18 348 

19 357 

20 369 

AVERAGE 336 



390 


-5.1 


330 


1. 8 


430 


-19. 1 


390 


- 8.5 


380 


- 2.9 


3 84 


-12.5 



indicates positive heat exchange 



prevented a drastic decrease in the forecasted layer depth thus averting 
the large errors which occurred at Station n N M as previously noted. 

A stated desire of this phase was to have a minimal difference 
between the final forecast layer depth and the observed condition. 

This was accomplished for this forecast period with a difference of 
only 11 feet (-2. 9%) on the last day. 

The average standard deviation for this period was 48 feet 
with an error of -12. 5% in forecasting the average layer depth. 

3. September, Station n N M (00Z forecasts) 

Results are shown in Table VIII. The results were unsatis- 
factory with errors ranging from 50% to 85% between the forecasted 
and observed layer depths. This was due again to light mechanical 
mixing which allowed the initial layer depth on the first day to decrease 
rapidly (140 feet to 29 feet) while the layer actually deepened 20 feet. 
This period was dominated by positive heat exchange (3 of 4 days) with 
convective mixing on the last day being ineffective in producing an 
accurate forecast due to the shallow initial layer depth. 

Again, the effect of forecasting persistence for the first 
three days with the calculation of the layer depth on the fourth day 
based on the persisted layer depth was investigated. The errors in 
this case between the forecast and observed layer depths were in the 
range, 13% to 21%. In predicting the average layer depth, this mod- 
ified model had an error of 2% and an average standard deviation of 
33 feet. 



49 



TABLE VIII. Comparison of the forecast MLD to the observed MLD 
(in feet) and the percent error between the two for 
Station M N n , September (00Z forecasts). 



DATE 


FORECAST 


OBSERVED 


% ERROR 


*11 


29 


160 


-81. 9 


*12 


19 


120 


-84. 2 


*13 


53 


no 


-51. 8 


14 


56 


180 


-68. 9 


AVERAGE 


39 


143 


-72. 7 



indicates positive heat exchange 



4. September, Station n N ,! (12Z forecasts) 



Results are shown in Table IX. This period also produced 
bad results with the exception of the forecast for the 13th which had a 
-1.7% error. Errors on the remaining days ranged from 45% to 70%. 

The first three days were controlled by the positive heat 
exchange and thus, with no convective mixing occurring, the layer 
depth was determined by mechanical mixing. Since the mixing was 

small (3 to 5 foot combined wave height) the surface heating produced 

✓ 

a shallow transient layer of 60 feet after three days and the only good 
forecast. When convective mixing occurred the following day, the 
transient layer was destroyed but the convective mixed layer depth 
model in use is such that a drastic increase in layer depth cannot be 
forecasted. Therefore, although convective mixing was the driving 
force during the final three days, the forecasted layer depth continued 
to lag behind observed conditions. 

5 . Summary 

Large deviations of the model from observed conditions 
occurred mainly when an amount of heat was added and mechanical 
mixing followed for one or more days causing the formation of a shal- 
low mixed layer (as shown by figures 3 through 5 of Appendix G); the 
convective mixing model was unable to quickly eliminate the shallow 
thermocline upon the return to a negative heat exchange. 

The one series of good forecasts occurred for a period when 
both convective and mechanical mixing were substantial. The large 
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TABLE IX. Comparison of the forecast MLD to the observed MLD 
(in feet) and the percent error between the two for 
Station "N", September (12Z forecasts). 



DATE 


FORECAST 


OBSERVED 


% ERROR 


*11 


39 


130 


-70. 0 


* 12 


49 


160 


-69. 3 


*13 


59 


60 


- 1.7 


14 


59 


160 


-63. 1 


15 


72 


130 


-44. 6 


16 


74 


180 


-58. 9 


AVERAGE 


59 


137 


-56. 9 



* indicates positive heat exchange 
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mechanical mixing prevented the formation of a shallow layer until 
convective mixing again became dominant. 

As in the daily forecasts discussed in the previous section, 
the need is evident for some method of predicting the formation of a 
shallow layer during periods when mechanical mixing is dominant. 
Possible methods of accomplishing this are discussed in the summary 
of the previous section. 

D. CONVERGENCE/ DIVERGENCE CORRECTIONS 

To determine the effectiveness of the convergence/divergence 
corrections, the daily forecasts were studied to determine whether or 
not the correction improved the accuracy of the forecasted layer 
depth. The daily forecasts were selected in order to avoid errors 
caused by poor input data. 

In determining whether the convergence/divergence corrections 
were of any advantage, any correction less than one foot was deter- 
mined to have produced no effect. Also, no determination of the effect 
of the correction on the accuracy of the forecast was made for those 
cases in which there was a large discrepancy between the forecast 
and observed layer depths. 

The first forecasts studied were for Station M N n in September and 
January. The maximum correction for this set of data was three feet. 
Of 15 forecasts which were considered to exhibit satisfactory accuracy, 
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convergence/divergence had no effect in nine cases. In the remaining 
six cases, the correction improved the accuracy half the time. 

The second forecast period studied was Station n P M in January. 
Four of the six corrections for this period were in the range, 9 to 1 1 
feet. Of the six forecasts, one correction produced no effect while 
three decreased the accuracy and two improved it. 

The final period investigated was January at Station !, E n . Four 
of the six corrections in this case were in the range, 11 to 19 feet. 
Applying the correction increased the accuracy in one case while it 
was decreased in two others. There were three forecasts in which 
the correction had no effect. 

It is evident that although there are locations where convergence/ 
divergence has very little effect (Station M N n ) there are areas where 
it must be taken into account (Stations n P n and n E n ). It was also noted 
that of 14 occasions when the convergence/divergence correction had 
some effect, it improved the accuracy of the forecast only in six 
cases. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



The results of this study indicate that the model which was 
developed will produce a slightly more accurate forecast than the 
FNW C formula when there is substantial convective action and the 
conditions vary from climatological conditions (Table II). Under 
these conditions the model had a maximum error of 7. 9% and a stan- 
dard deviation of 29 feet. This is compared to a maximum error of 
12. 5% and a standard deviation of 53 feet for the FNW C formula. 

The greatest source of error in the model developed in this study 
is when convective mixing ceases for a short time. When this happens, 
the forecast layer depth is based on the layer depth due to mechanical 
mixing and this at times forecasts a decrease in the layer depth of up 
to 400 feet. This is a serious flaw in the model and work needs to be 
done in this area. What is actually happening, is the formation of a 
new shallow layer (rather than the deep layer being destroyed) for a 
short period of time with a return to a deep layer as soon as either 
convective or mechanical mixing becomes sufficiently strong to destroy 
this transient shallow layer. Work must be done in this area to deter- 
mine under exactly what conditions a shallow layer will form or be 
prevented from forming. Two methods which might be used to improve 
the accuracy of the forecast until such time as improved methods are 
available are as follows: 
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1. Make a persistence forecast in the event of a forecast of a 
drastic change in layer depth. 

2. Allow the layer depth to decrease only a certain percent of 
the difference between the previous layer depth and the fore- 
cast layer depth. 

In comparing the mechanical mixing relationships, the Laevastu 
method (used in this study) was more accurate than the FNWC method. 

In the one good forecast of mechanical mixing investigated, the FNWC 
relation had an error of 25% (compared to 1. 7% for the Laevastu 
relation) and the next best result obtained was an error of 29 . 3% 
(compared to 18.2% for Laevastu's method). An area of further study 
concerns this Laevastu method. It is possible that the constant should 
be changed, perhaps to a second order term in combined wave height. 
This is suggested by the fact that as the combined wave height increased 
(from 5 to 27 feet) the error between the forecast and observed con- 
ditions decreased (from 70% to 18%). 

It appears that convergence/divergence is a factor, particularly 
for Stations "P" and n E n . However, the method used in this study to 
compute the correction for convergence/divergence produced a correc- 
tion which improved the accuracy of the forecast less than half the 
time. Since FNWC is now in the process of improving their surface 
current analysis model, this factor should be investigated at a future 
time. 
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The basic conclusion reached from this study is that the James 
method for convective mixing can be computerized and forecasts can 
be produced which take the heat budget into consideration. The ac- 
curacy of this model is restricted by the accuracy of the input data, 
particularly that of the initial layer depth. Given an accurate layer 
depth, this model has demonstrated the ability to make a continuing 
forecast using forecast layer depths as input data (during periods of 
large negative heat exchange) until such a time as a new BT report is 
available (Table VIII). 
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APPENDIX A 



FNWC POTENTIAL LAYER DEPTH METHOD 

The first step in the FNWC layer depth forecast method is the 
construction of the layer depth analysis. 

The initial step in the analysis is to adjust the 12-hour-old 
potential layer depth field toward the interpolated climatological field 
of the analysis time. This interpolated field is a linear interpolation 
over two months; the present and preceding month if the day is less 
than 15, or the present and following month if the day is more than 15. 
If the analysis is for the 15th, only the climatology of the present m 
month is used. This adjusted climatological field is the first guess 
field. 

Next, the mixed layer depth possible due to wave mixing alone is 
computed from the combined wave height analysis as follows: 

MLD = CH[ 11-0. 1(SST-T 6oo )] sst * t 6oo 

= 11 CH , SSTc T . 

600 

The effect of the term, 0. I(SST-T^qq), is to reduce the depth of 
mixing when a strong thermocline gradient exists. ‘ 

The mixed layer depth field due to wave action is then smoothed 
with a special smoothing function designed to decrease the highest 
values in the field without increasing the values at any point. 
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Next, the MLD field due to wave mixing is compared to the 
initial-guess field produced from climatology and the deeper of the 
two is chosen as the second-guess MLD field. 

The final step in the analysis phase is to adjust this second guess 
field to BT reports taken within the past 72 hours. The reports are 
weighted, based on the age of the observations at analysis time with 
the latest reports receiving the greatest weight. One final control 
placed on the analysis field is to force the layer depth to be less than 
or equal to 900 feet at all grid points. The result is the FNWC poten- 
tial mixed layer depth analysis field. 

The forecast layer depth field is produced in the same manner as 
the analyzed field. The analyzed field is adjusted toward the clima- 
tology of the forecast time using the method described for producing 
the analysis. The mechanical mixing is then computed using the 24- 
hour forecast combined wave field. The same steps of comparison and 
smoothing are taken between the wind mixed field and the clima- 
tological field. The result is the forecast potential mixed layer depth 
field. 

The one difference in procedure between the analysis and fore- 
cast phase is that the forecast field is not adjusted to BT reports. 

A severe constraint on this method is the lack of data which 
results in relying essentially on climatology in areas of sparse data. 

An additional problem area is that this forecast scheme is also subject 
to errors in the wave height forecast. 
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APPENDIX B 



FNWC SEA SURFACE TEMPERATURE ANALYSIS 

The sea surface temperature analysis is accomplished through a 
scheme called Fields By Information Blending (FIB). The analysis 
is done on a 125x125 grid with an output on a 63x63 grid. 

Inputs to the analysis are sea surface temperature reports within 
six hours of the analysis time and surface temperatures from BT 
reports received in the past 12 hours. When it becomes readily 
available, satellite data also will be incorporated in the model. 

The first step in the analysis is to compute the current clima- 
tology. This climatology is then assigned a weight based on the 
gradient of temperature between grid points (highest gradients are 
given the lowest weights). Next, the first guess field is produced by 
blending the previous analysis and the current climatology. 

The next step is to read in and screen the sea surface tem- 
perature reports. A difference field is computed by subtracting the 
first guess field from the field produced from reports. Each report 
is assigned a weight based on its age. 

The final sea surface temperature analysis is produced by blend- 
ing this difference field with the guess field using the assigned weights. 
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The basic drawback to this procedure is the fact that 1100 to 1200 
reports are input to produce an analysis over 3969 grid points. Ob- 
viously, some accuracy is lost through smoothing. The procedure 
should become better once satellite data are routinely available for 
use as input parameters. 
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APPENDIX C 



FNWC SUBSURFACE THERMAL STRUCTURE ANALYSIS 

The subsurface temperature analysis is used for the following 
standard levels: surface, 100 feet, 200 feet, 300 feet, 400 feet, 600 

feet, 800 feet, and 1200 feet. 

The inputs to the model are: the 12-hour-old temperature analysis 

at standard levels, climatology at standard levels, sea surface tem- 
perature analysis, mixed layer depth analysis, and BT observations 
for the past 72 hours. 

The first step is to prepare guess fields for each standard level. 
This is done in two steps. First the 12_hour-old analysis is adjusted 
toward climatology; then this field is adjusted to the sea surface tem- 
perature and mixed layer depth analyses to produce the final guess 
field. 

The second step is to prepare data lists from reported observa- 
tions. Each report is assigned a weight factor on the basis of age and 
then interpolated to find temperatures at all standard levels from the 
surface to the bottom of the report. 

The final step is to adjust the fields to the data. First the final 
guess field is adjusted to the data and a light smoothing function is 
applied to the vertical gradient between the level in question and the 
level above. This results in a first analysis field. Next, the first 
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analysis field is adjusted to the data and a vertical control is applied 

requiring that the temperature at the level in question is less than or 

equal to the temperature of the level above. The final analysis for a 

level is produced after the temperature at any level is restricted to 
o o 

the range -2 C to 35 C. 

The primary limitation in this method is the lack of data. In 
addition, areas of strong currents and areas of strong winds which 
are not forecast will produce errors in the forecast. 
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APPENDIX D 



FNWC HEAT EXCHANGE MODEL 

The total heat exchange across the air- sea interface is a com- 
bination of five components: heat absorbed from the sun, solar heat 

reflected from the sea surface, longwave radiation from the sea sur- 
face less the incoming longwave radiation from the atmosphere, 
sensible heat exchange across the air-sea interface, and latent heat 
exchange thru evaporation or condensation. Each component was 
calculated from a separate equation and the results were combined to 
obtain the net heat exchange. Following are the equations used in the 
computation of each component: 

1. Solar insolation, Q 

s 

3 

Q = 0. 014 A tJl. 0 - 0. 0006 C ) 
s n d 

2. Reflected solar insolation, Q 

r 

Q = 0. 15Q - (0. 01Q ) 2 

r s s 

3. Effective back radiation, Q, 

b 

Q b = Q b (l. 0 - 0. 765 C) where, 

Q ob = <rT l [A+4( . T w - Ta f] where, 

3 -4 

A = 0. 261 exp [-7. 77 x 10 (273-T )] 

a 

4. Latent heat transfer, Q 

e 

Q = (0.26+0. 077V) (0.98e-e a )L f [ e -e_ positive] 

Q e = 0. 077V(0. 98e w _ e a )Lt t e w- e a negative] 
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5. Sensible heat transfer, A 



Cl = 39(0. 26 + 0. 077) (T _T ) 
h w a 

Q h = 3V(T„ - T a > 



[T -T positive] 

w a r 

[T -T negative] 

L w a 



These five parameters are then combined to give total heat exchange, 
Q , as follows: 



These formulas were selected through a comprehensive literature 



of performance in studies supported by field testing. 

The accuracy of these formulas in synoptic use can be determined 
directly only in rare instances when a research vessel is available to 
make radiation measurements. However, it is possible to make an 
indirect verification by comparing the heat exchange analyses with 
actual measured changes which occur in the ocean and atmosphere. 
This type of indirect verification has been used for the past several 
years, and there has been no evidence that the FNW C heat exchange 
formulas are inaccurate or in need of change. 

The errors which do occur in the heat exchange computations are 
probably caused by inaccuracies in the input analyses and forecasts. 

In general, these input parameters are reasonably accurate but errors 
in meteorological observations and in the manipulating and averaging 
of sparse data over large sea areas can cause problems. Cloud 
analyses have been one of the greatest error sources. The wind field 
analyses are generally good, but a small error in the wind field at 



n 




review and were those which promised the greatest reliability because 



65 



high wind speeds will cause large errors in some heat exchange 



components, particularly in latent heat exchange. 



66 



APPENDIX E 



FNWC SURFACE CURRENT ANALYSIS 

Ocean currents are caused by the combined action of several 
physical environmental forces. Thermohaline circulation, wind and 
mass transport by waves, and tidal effects are the most important 
of these forces. Other forces such as changes in atmospheric pressure, 
inertia currents, current boundaries, and the coriolis effect were 
considered, but were discarded for lack of significance or for lack of 
any forecasting method. Tidal currents are important in coastal 
waters, and have some magnitude in deep water. However, it was 
decided that for a 24-hour prediction, the tidal currents could be 
neglected. The elimination process finally led to the assumption for 
the model that the total current flow could be described by the vector 
sum of permanent thermohaline flow and local wind-driven flow, 
including mass transport by waves. 

The thermohaline or permanent circulation evolves from the 
distribution of ocean temperature and salinity. The thermohaline- 
type currents have been studied for some 70 years, and have been 
modeled by geostrophic methods. A drawback in applying this geo- 
strophic method is that it requires accurate information on the dis- 
tribution of temperature and salinity in the water column for which 
little real-time salinity data is available. Therefore, the direct 
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application of the geostrophic method at FNWC has been modified to 
calculation of a permanent flow component based on the local tem- 
perature structure alone. 

The permanent component calculation starts by computing a 
weighted average temperature from the surface to the 600-foot depth. 
Values for the two levels are taken from FNWC synoptic analyses. 
Using this average temperature, the thermal wind equation is applied 
to obtain u and v components of the surface current caused by distribu- 
tion of sea water temperatures. 

The second component in the surface current model is the local 
wind-driven flow. Complicated calculation methods are available, 
but they are not well adapted to synoptic analysis so an empirical 
approach has been used. A simple formula (a constant times the wind 
speed) was developed relating wind to surface currents. The constant 
is adjusted to account for mass transport by waves. The wind is 
computed from the FNWC sea level pressure analysis and the 24-hour 
sea level pressure forecast. The geostrophic wind is corrected for 
stability and curvature to obtain a "marine" wind. 

Thus, the FNWC total ocean surface current is the vector sum of 
the thermal and wind-driven components. 

The surface current model is designed to produce an accurate 
picture of current patterns over deep water. It will do this, but it 
makes no pretense of accuracy in other areas and for all possible 
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situations. The grid size of about 200 miles precludes accuracy in 
areas of marked thermal gradients, such as along the edge of strong 
currents like the Gulf Stream. Very narrow currents can pass between 
grid points without showing in the surface current analysis. The model 
does not perform well where for those few currents which are dom- 
inated by a haline distribution, as in areas where there is much melt- 
ing of ice or alternate periods of melting and freezing. A final area 
of model weakness is in coastal areas or near island chains. 
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APPENDIX F 



FNWC SINGULAR ADVECTIVE WAVE AND SWELL MODEL 

The total n sea state" may be described as a spectrum of sinusoidal 
waves of varying frequency and direction of movement. 

The singular model typically relates wave height to wind speed 
through a simple exponential r elationship. It is significant to note 
that singular models usually do not maintain wave continuity from one 
computational cycle to the next, nor is there provision for energy 
sinks through dissipation due to opposing winds. 

Spectral models treat the growth, propagation and dissipation of 
energy for a number of frequencies traveling in a number of directions. 
These models divide the frequency and direction spectrum into a finite 
set of bands, calculate the changes in energy for each frequency- 
direction component and integrate over-all components to obtain energy 
moments which statistically describe the total sea state. More 
sophisticated spectral models contain wave/wave interaction terms in 
addition to those relating wind/wave interaction. 

Since 180 components are treated at each grid point, spectral 
models are very demanding on computer time. The singular advective 
model presently in use represents a compromise between computer 
limitations and model completeness. It utilizes the propagation of 
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energy concept but only the frequency of maximum energy is propogated 
in the direction of the wind. (A single component is selected as being 
representative of the total spectrum. ) Growth and dissipation are 
achieved through comparison of advective heights with fully developed 
heights from a new wind field at a later time. 

The starting points for each analysis cycle are the six-hour old 
wind wave and swell fields. The previous wind wave analysis is 
advected forward in two three-hourly time steps with dissipation and/ 
or growth through comparison with the current wind analysis. The 
product of these two advective steps is a preliminary calculation of the 
new wind wave analysis. 

A final wind wave analysis is obtained through the introduction of 
synoptic ship reports into the preliminary calculation of wave heights. 

The forecast portion of the model is simply a continuation of the 
advective propagation outlined above except that forecast winds are 
used in the dissipation and growth of the wind waves. 

The six-hour-old swell trains are similarly advected forward, 
but without consideration of any interaction with current winds. 

(This is a weakness which is not inherent to true spectral models 
where wind/wave and wave/wave interactions occur continually for 
all frequencies. ) 

The speed of propagation is proportional to swell period. Since 
this model makes no provision for the dissipation of swell energy 
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through swell/wind interaction, swell heights could propagate unchanged. 
To insure that there will always be some decay in swell heights with 
time, a decay coefficient has been added to the advective equations. 

There are three types of limitations which contribute to inac- 
curacies in sea/swell analyses and forecasts: 

a. Wind errors. 

b. Wave observation errors. 

c. Modeling limitations. 

Even if model formulations were perfect, wind errors would lead to 
incorrect wave computations. Determination of the wind at the surface 
of the water which actually controls wave growth is particularly 
difficult. The forecast winds used in the forecast cycle also contain 
any errors resulting from deficiencies in the model used to produce 
the winds. 

Wave observation errors also contribute to errors in the final 
wind wave analysis. This is a significant factor since almost all 
wave height reports represent an ‘'eyeball” estimate. 

The inability of any singular model to define the complete wave 
energy spectrum and all its interactions is a major modeling limitation. 
The advective singular model represents a significant improvement 
over earlier singular approaches for it insures wind wave continuity 
between runs and treats energy dissipation due to opposing winds. 
However, wave/wave interactions which influence both growth and 
dissipation are not considered at all. 
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APPENDIX G 



SELECTED BT TRACES 

Figures 3 thru 5 show another example of the formation of a 
shallow surface layer. The difference between this case and the 
preceding one is that in this example the layer is formed not due to 
heating but due to wind mixing during a period of slight negative heat 
exchange. Figure 5 shows the return to a deep layer following & 
period of substantial negative heat exchange. 



73 



DEPTH (feet) 



TEMPERATURE (°F) 




Figure 3. BT trace for 00Z 18 September, Station "N". 
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Figure 4. BT trace for 00Z 19 September, Station "N". 
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Figure 5. BT trace for 00Z 20 September, Station "N". 
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APPENDIX H 



PROGRAM FOR DAILY MIXED LAYER DEPTH FORECASTS 



DIMENSION I T ( 4) , I TS (4 ) , I QN ( 4 ) , I U 1 4 ) , IV (4 ) , I WH (4 ) , 

ISH(4) ,U(4),V (4 ) 

N = 0 

WRITE! 6,5 01) 

501 FORMAT! • 1», ////////, 27X, 'FORECAST MLD FOR STATION PAPA*) 
9000 N=N+1 

IF(N.EQ.9) GO TO 9999 
IFIN.EQ.6) GO TO 5000 
GC TO 5 5C0 
500C WRITE!6,500L) 

5001 FORMAT ( ' 1 ' ,// //////, 27 X, 'FORECAST MLD FOR STATION PAPA 
{ CCNT ) ' ) 

55CC READ! 5,500) KT.KD 

R EAD!5 , 1 ) ( IT ( 1 ) , 1 = 1 ,4) 

READ! 5,1) ( ITS! I ) ,1 =1 ,4) 

READ! 5, 1) ! I QN ( I ) , 1=1, 4) 

R EAD! 5 , 1 ) ! 1U ( I ) , 1=1 ,4 ) 

READ! 5,1) ( I V! I ) , I =1,4) 

R EAD! 5, 1 ) ( IWH( I ) , 1=1,4) 

R E AD! 5 , 1 ) (I SH( I ) , 1=1 ,4) 

R EAD! 5,2) PL I 
READ! 5 , 2 ) RLD 

1 FORMAT (416) 

2 FORMAT ( F 7 . 2 ) 

500 FORMAT (214) 

INTERPOLATION FACTORS X C Y 
X = 

Y = 

INTERPOLATE TEMPERATURE 

T L= IT ( 1)-X*( IT! 1 ) - I T { 2 ) ) 

TU = I T { 3 ) — X" ( IT (3 )— -IT{4) ) 

TF = TU- Y ; '( TU-TL) 

INTERPOLATE TEMPERATURE 600 

TSL=ITS(1)-X~(ITS(1)-ITS(2) ) 

T SU= I T S ( 3 J — -X ~'( ITS! 3)— ITS! 4) ) 

T SF=T SU-Y : * ( TS U-T S L ) 

INTERPOLATE HEAT EXCHANGE 

. QNL= I QN ( L ) -X" ( IQM ( l)-IQN(2) ) 

QNU = I QN (3 ) -X* ( I QN (3 ) - I QN (4 ) ) 

QNF=QNU-Y* ( QNU-QNL ) 

INTERPOLATE WAVE HEIGHT 

WH L = I WH ( 1 ) -X” (IWH(l) — IWHC2 ) ) 

W HU= I WH ( 3 ) - X-- ( IWH< 3 ) - I WH(4) ) 

W F F = WHU-Y * ( W HU-W h L ) 

INTERPOLATE SWELL HEIGHT 

SHL= I S H ( l)-X*(ISH{ I )— I SH ( 2 ) ) 

SHU= I SH (3 )-X"(ISH(3)-ISH(4) ) 

SFF=SHU-Y J ~(SHU-SHL) 

REDUCE DATA BY SCALING FACTOR 
TC=TF/64. 

TSC=TSF/64. 

. T={ (9.*TC)+160. >/5. 

TS=((9.*7SC)+160.)/5. 

QN=QN F/ 16 . 

Wh=WHF/512 . 

SH=SHF/ 512 . 

DO 100 J = 1,4 
U(J)=IU(J)/16. 

V ( J ) = I V ( J ) / 16 . 
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100 CONTINUE 
CP=0. 935 
RW= 1.024 
DCCN=0 .0 

IF(QN.GE.O.O) GO TO 1200 
CONVECTIVE MIXING 

DT = ( 1 00 .* ( T-T S ) ) / ( PL 1-600 . ) 

CM=! 2. '-CN ) /( CP*RW'*DT) 

CCCN= (PLI**2+CM) *-0 . 5 
MECHANICAL MI XING 
1200 CH=!WH**2+SH**2)**0.5 
DM = 12 .5*CH 

I F(DM.GE.OCON) GO TO 1300 
DL = DCQN 
GO TO 1000 
1300 DL=OM 
CON V/D I V CORRECTION 

1000 A=U(1 ) +U ( 3 ) — U ( 2 ) — U (4 ) +V ( 3 ) +V (4 ) —V ( 1 )-V ( 2 i 
GRID SPACING *G' 

G= 

B = ( A-~0 .864 )/ !4 . *G ) 

DO=B"D L 
PRLD= DL +DD 
WRITE(6,502) KT,KD 
WRITE! 6,3000) DCON 
WRITE (6, 3 100) DM 
WRITE(6 ,3200 ) DL 
WRITE! 6 ,3300) DD 
WRITE (6, 3400 ) PRLD 
KPI TE ( 6 ,3 5 00) RLD 

502 FORMAT!' • , /// / , 1 3 X , • F OREC A ST FOR ',12,'Z ',12,' 
3000 FORMAT! 'O' ,23X, 'CONVECTIVE MLD IS ',F8.2) 

3 1 CC FORMAT!' ' ,23 X , • F CRC E C MLD IS ',F8.2) 

3200 FORMAT!' ',23X,' LAYER DEPTH IS *,F8.2) 

3300 FORMAT ( * ' , 23X , ' CONV/ 0 I V CGRR IS ',F8.2) 

3400 FORMAT!' • ,23X , • F ORECAST MLD IS ',F8.2) 

3500 FORMAT!' ', 23X ,' ACTUAL MLD IS ',F8.2) 

GO TO 9000 
9995 STOP 
ENC 



JAN.') 
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PRCGRAM FOR CONTINUOUS MIXED LAYER DEPTH FORECASTS 



DIMENSION IT (4) , ITS (4 ) , IQN(4 ) , IU (4 ) , IV ( 4 ) , IWH(4), 

I SHU) , U(4) ,V(4) 

N = 0 

InPITEU ,501 ) 

501 FORMAT! » 1 ',////////, 22X , 'CONTINUOUS MLD FORECAST FCR 
STATION PAPA' ) 

R t AD ( 5 » 2 ) PLI 
PRLD=Pl I 
9000 N = K + 1 

IFIN.EQ.9) GO TO 9999 
IFIN.EQ.6) GO TO 5000 
GO TO 5500 
50CC Vs R I TE ( 6,5001) 

5001 FORMAT { * l* ,////////, 27X, 'FORECAST MLD FOR STATION PAPA 
( C 0 NT ) ' ) 

5 5 OC READ! 5, 500 KT , KD 

R E AD ( 5 , l) ( IT { I) , 1=1,4) 

R E AD ( 5 , 1 ) (ITS( I) ,I=L ,4) 

R E AC( 5,1) (IQN( I ) , 1 = 1,4) 

R E AD ( 5 , I ) (IUt I), I = L,4) 

READ (5,1) ( I V ( I ) ,1=1,4) 

READ! 5, l) ( 1 W H ( I ) , I = L , 4 ) 

R E A C ( 5 , 1 ) ( ISH( I ) , 1 = 1,4) 

PLI =PRLD 
R E A D ( 5 , 2 ) PLD 

1 FORMAT (416) 

2 FORMAT (F7.2) 

500 FORMAT (2145 

INTERPOLATION FACTORS X & Y 
X = 

Y = 

INTERPOLATE TEMPERATURE 

T L= I T ( 1 )-X*( I T ( 1 ) — I T ( 2 ) ) 

TU= IT (3 )-X*( IT (3 )- IT ( 4 ) ) 

TF=TU-Y*( 1 U-TL) 

INTERPOLATE TEMPERATURE 600 

TSL=ITS(1)~X*(ITS(1)-ITS(2) ) 
TSL=ITS(3)-X«(ITS(3J-ITS(4) ) 

T S F=T SU-Y = ( T SU-T SL ) 

INTERPOLATE HEAT EXCHANGE 

QNL=I QN ( i)-X*(IQN(U -I ON ( 2 ) ) 

•• CNU=1 QN ( 3 ) -X* ( I QN (2)-IQN(4) ) 
t NF=QNU-Y ;< ( CNU-QNL ) 

INTERPOLATE WAVE HEIGHT 

W h L = I W H ( l)-X*( IWH ( l )- IWH( 2 ) ) 

WHU=I WH ( 3 ) -X* (IWH(3)-lWh(4) ) 

NFF=WHU-Y*( WHU-WHL ) 

INTERPOLATE SWELL HEIGHT 

SHL=ISH(1 ) -X* ( I SH ( 1 )-ISH(2) ) 

SHU=ISH(3)-X*( ISH(3)-ISH(4) ) 

5 FF=S HU-Y ~ ( SHU-S H L ) 

REDUCE DATA 3Y SCALING FACTOR 
TC=TF/64 . 

T SC=TSF/64 . 

T = ( ( 9 . * T C ) + 1 6 0 . ) / 5 » 

TS=((9.*TSC)+160.)/5. 

CN=CNF/16. 

Wh= WHF / 512. 

Sh=SHF/ 512 . 

DC 100 J = 1 ,4 
U ( J )= I U ( J ) / 16. 

V( J ) = I V ( J ) / 16 . 

ICC CONTINUE 
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CP = 0 .935 
R Vs = 1 . 02 4 
DC0N= 0 . 0 

IF(CN.GE.O.O) GO TO L200 
CONVECTIVE MIXING 

DT=( LOO.*(T-TS) )/( PL 1-600. ) 

C M = ( 2 .*CN)/(CP-'RW*DT ) 

DCQN=( PLI**2+CM) **0. 5 
MECHANICAL MIXING 
12 CO CH=( WH**2 + SH**2 )**0.5 
DK=12.5*CH 

IF(CM.GE.DCON) GO TO 1300 
DL=DCON 
GO TO 1C00 
L300 CL=CM 
CCNV/DIV CORRECTION 

IOOC A=U( I) + U(3)-U(2)-U(4)+V(3) + V{4)-V( 1 ) — V{ 2) 

GRID SPACING • G' 

G = 

B=(A*0.864)/(4.*G) 

DC=B*DL 
PRLO=DL+DD 
WR I TE ( 6 » 502) KT , KD 
WR IT E ( 6 » 3000 ) DCON 
WRITE(6,3L00) DM 
WRITE ( 6,32001 DL 
WR ITE (6 ,3300 ) DO 
WRITE! 6 ,3400) PR LD 
WR I TE ( 6 , 3500 I RLD 

502 FORMAT ( * ',////, L3X FORECAST FOR *,I2,*Z ',12, 

30CC FORMAT! '0* ,23X, 'CONVECTIVE MLD IS ',F8.2) 

3 LOO FORMAT!* ' , 2 3X , • Fu RC E C ML D IS «,F8.2) 

32 OC FORMAT ( • *,23X,' LAYER DEPTH IS *,F8.2) 

330C FORMAT!' ' , 2 3X , • CON V/D I V CORR IS *,F8.2J 
3400 FORMAT!' • , 23X , « FOR ECAST MLD IS *,F8.2) 

3 5 CC FORMAT! » ' ,23 X , • AC TUAL MLD IS »,F8.2) 

GO TO 9000 
9999 STCP 
END 



JAN . •) 
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APPENDIX I 



FORECAST MLD FOR STATION NOVEMBER 



FGRECAST FOR OZ 13 SEP. 

CONVECTIVE MLD IS 0.0 

FORCED ML C IS 53.71 
LAYER DEPTH IS 53 .71 
CONV/DIV CORR IS -0.37 

FORECAST MLD IS 53.34 
ACTUAL MLD IS 110.00 



FORECAST FOR 12Z 13 SEP . 

CONVECTIVE MLD IS 0 .0 

FORCED MLD I S 59.25 

LAYER DEPTH IS 59.25 
CONV/DIV CORR IS -0.11 

FORECAST MLD IS 59.14 
ACTUAL MLD IS 60.00 



FORECAST FOR OZ 14 SEP. 

CONVECTIVE MLD IS 110.26 
FORCED MLD IS 55.67 

LAYER DEPTH IS 110 .26 

CONV/DIV CORR IS -0.17 
FORECAST MLD IS 110.09 
ACTUAL MLD IS 180.00 



FORECAST FOR 12Z 14 SEP. 

CONVECTIVE MLD IS 60.21 
FORCED MLD I S 55. 18 
LAYER DEPTH IS 60. 2L 
- CONV/DIV CORR IS -0 . 1 1 
FORECAST MLD IS 60.10 
ACTUAL MLD IS L60.00 



FORECAST FOR 12Z 15 SEP. 

CONVECTIVE MLD IS 160.49 

FORCED MLD IS 73.91 

LAYER DEPTH IS 160.49 

CONV/DIV CORR IS 1.09 

FORECAST MLD IS 161.58 
ACTUAL MLD I S 130. 00 
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FORECAST MID FOR STATION NOV EM BE R ( CONT ) 



FORECAST FOR OZ 16 SEP. 

CONVECTIVE MLD IS 150.24 

FORCED MLD IS 51.32 

LAYER DEPTH IS 150.24 

CONV/DIV CORR IS 1.11 

FORECAST MLD IS L51.35 
ACTUAL MLD IS 150. OC 



FORECAST FOR 12Z 16 SEP. 

CONVECTIVE MLD IS 130.50 
FORCED MLD IS 53.70 

LAYER DEPTH IS 130.50 

CONV/DIV CORR IS 1. 01 

FORECAST MLD IS 131.51 
ACTUAL MLD IS 180.00 



FORECAST FOR OZ 15 SEP. 

CONVECTIVE MLD IS 190.11 

FORCED MLD IS 37.08 

LAYER DEPTH IS 190.11 

CONV/DIV CORR IS 1.03 

FORECAST MLD IS L91.14 
ACTUAL MLD IS 20.00 



FORECAST FOR 12Z 19 SEP. 

CONVECTIVE MLD IS 160.08 
FORCED MLD IS 32.67 

LAY ER DEPTH IS 160 .08 

. CONV/DIV CORR IS 0. 11 

FORECAST MLD IS 160.19 
ACTUAL MLD IS 140.00 



FORECAST FOR OZ 20 SEP. 

CCNVECT IV E MLD IS 22 .52 

FORCED ML C IS 27.62 
LAYER DEPTH IS 27.62 
CONV/DIV CORR IS 0.21 

FORECAST MLD IS 27.83 
ACTUAL MLD IS 170.00 
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FORECAST MLC FOR STATION NOVEMBER (CCNT ) 



FORECAST FOR L2Z 2C SEP. 

CONVECTIVE MLD IS 140.22 

FORCED MLD IS 29.36 

LAVER DEPTH IS 140.22 

CCNV/DIV CORR IS 0.50 

FORECAST MLD IS 140.72 
ACTUAL MLD IS 170.00 
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FORECAST MLD FOR STATI CN NOVEMBER 



FORECAST FOR OZ L2 JAN. 

CONVECTIVE MLD IS 550.06 
FORCED MLD IS 132.37 
LAYER DEPTH IS 550.06 
CONV/DI V CORR IS I .90 

FORECAST MLD IS 551.96 
ACTLAL MLD IS 530.00 



FORECAST FOR L2Z 12 JAN. 

CONVECTIVE MLD IS 560. 04 
FORCED MLD IS 125.29 
LAYER DEPTH I S 560.04 
CONV/DIV CORR IS 1.44 

FORECAST MLD IS 561.48 
ACTUAL MLD IS 570. CO 



FORECAST FCR OZ 13 JAN. 

CONVECTIVE MLD IS 530.04 
FORCED ML D I S 138.08 
LAYER DEPTH IS 530.04 
CONV/DIV CORR IS 0.76 

FORECAST MLD IS 530.80 
ACTUAL MLD IS 510.00 



FORECAST FOR 12Z 13 JAN. 

CONVECTIVE MLD IS 570.02 
FORCED ML C IS 73.65 

LAYER DEPTH IS 570.02 

CCNV/DIV CORR IS 1.94 

FORECAST MLD I S 571 .95 
ACTUAL MLD IS 530. CO 
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FORECAST MLD FOR STATION PAPA 



FORECAST FOR OZ 13 JAN. 

CONVECT IV E MLD IS 390. 15 
FORCED ML C IS 164.92 
LAYER DEPTH IS 390. 15 
CONV/DIV CORR IS L0.L4 

FORECAST MLD IS 400.29 
ACTUAL MLD IS 360. CO 



FORECAST FOR 12Z 13 JAN. 

CONVECTIVE MLD IS 0.0 

FORCED MLD IS 130.68 
LAYER DEPTH IS 130 .68 
CONV/DI V CORR I S -0. 88 
FORECAST MLD IS 129.80 
ACTUAL MLD IS 430.00 



FORECAST FOR OZ 19 JAN. 

CONVECTIVE MLD IS 430.82 
FORCED ML C IS 139.04 
LAYER DEPTH IS 430.82 
COKV/DIV CORR IS 8.90 

FORECAST MLD IS 439.72 
ACTUAL MLD IS 390.00 



FORECAST FCR 12 Z 19 JAN. 

CONVECTIVE MLD IS 382.14 
FORCED MLD IS 171. 10 
LAYER DEPTH IS 382 .14 
•CONV/DIV CORR IS 3.53 

FORECAST MLD IS 385.67 
ACTUAL MLD IS 410.00 



FORECAST FOR OZ 20 JAN. 

CONVECTIVE MLD IS 391.68 
FORCED ML C IS 172.31 
LAYER DEPTH IS 391.68 
CONV/DIV CORR IS 10.94 
FORECAST MLD IS 402.63 
ACTUAL MLD IS 380. CO 
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FORECAST MLD FOR STATION PAPA(CCNT) 



FORECAST FOR L2Z 20 JAN. 

CONVECTIVE MID IS 410.00 
FORCED MLD I S L34. 26 
LAYER DEPTH IS 4L0.00 
CONV/DIV CORR IS 9.59 

FORECAST MLD IS 419.59 
ACTUAL MLD IS 460.00 
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FORECAST MLD FOR STATION ECHO 



FORECAST FOR OZ 12 JAN. 

CONVECTIVE MLD IS 858.16 
FORCED MLD I S 83.52 

LAYER DEPTH IS 8S8.16 

CONV/DIV CORR IS 11.99 
FORECAST MLD IS 910.15 
ACTUAL MLD IS 110.00 



FORECAST FOR 12Z 12 JAN. 

CONVECTIVE MLD IS 649.67 
FORCED MLD IS 92.04 

LAYER DEPTH IS 649.67 

CONV/DIV CORR IS 11. 13 
FORECAST MLD IS 660.81 
ACTUAL MLD IS 580.00 



FORECAST FOR OZ 15 JAN. 

CONVECTIVE MLD IS 769.89 
FORCED MLD IS 33.51 

LAYER DEPTH IS 765.89 

CONV/DIV CORR IS -0.14 
FORECAST MLD IS 769.75 
ACTUAL MLD IS 620.00 



FORECAST FOR 12Z 15 JAN. 

CONVECTIVE MLD IS 659.92 
FORCED ML C IS 52 .49 
LAYER DEPTH IS 659.92 
CONV/DIV CORR IS 18.80 
FORECAST MLD IS 678 .72 
ACTUAL MLD IS 650.00 



FORECAST FOR OZ 16 JAN. 

CONVECTIVE MLD IS 0 .0 

FORCED MLD I S 45. 71 
LAYER DEPTH IS 45.71 
CONV/DIV CORR IS 1.30 

FORECAST MLD IS 47.01 
ACTUAL MLD IS 100.00 
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FORECAST MLD FOR STATION ECHO(CCNT) 



FORECAST FOR L2Z 16 JAN. 

CONVECTIVE MLD IS 649.39 
FORCED MLD IS 68.38 

LAYER DEPTH IS 649.89 

CGNV/DIV CORR IS IS .86 
FORECAST MLD IS 665.75 
ACTUAL MLD IS 670.00 
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CONTINUOUS ML D FORECAST FOR STATION NOVEMBER 



FORECAST FOR 12 Z 1L SEP. 

CONVECTIVE MLO IS 0.0 

FORCED MLO IS 39.83 
LAYER DEPTH IS 39.83 
CONV/DIV CORR IS -0.73 

FORECAST MLO IS 39.10 
ACTUAL MLD IS 130.00 



FORECAST FOR L2Z 12 SEP. 

CONVECTIVE MLD IS 0.0 

FORCED MLD I S 49. 13 
LAY ER DEPTH IS 49 . 13 
CONV/DIV CORR IS -0.41 

FORECAST MLD IS 48.72 
ACTUAL MLD IS 160.00 



FORECAST FOR 12Z 13 SEP. 

CONVECTIVE MLD IS 0.0 

FORCED MLD IS 59 .25 

LAYER DEPTH IS 59.25 
CONV/DIV CORR IS -0.11 

FORECAST MLD IS 59.14 
ACTUAL MLD I S 60. 00 



FORECAST FOR 12Z 14 SEP. 

CONVECTIVE MLD IS 59.36 
FORCED MLD I S 55. 18 
LAYER DEPTH IS 59.36 
CONV/DIV CORR IS -0.11 
FORECAST MLD IS 59.25 
ACTUAL MLD IS 160.00 



FORECAST FOR 12Z 15 SEP. 

CONVECTIVE MLD IS 60.86 

FORCED MLD IS 73.91 
LAYER DEPTH IS 73. 9L 
CONV/DIV CORR IS 0.50 

FORECAST MLD IS 74.41 
ACTUAL MLD IS 130.00 



89 



Ur* 



FORECAST MLD FOR STATION NOV EM36 R ( CONT ) 



FORECAST FOR 12Z 16 SEP. 



CONVECTIVE MLD 
FORCEC MLD IS 
LAYER DEPTH IS 
CONV/DIV CORR IS 
FORECAST MLD IS 
ACTUAL MLD IS 



IS 75.38 
53.70 
75.33 
0.58 
75 .97 
180.00 



ic ^ r-s *i A ** A*: ^crcxc^ :& ^ ^3p £:£: si tH it *c •& 

^ ^ ^ ^ ^ ^ ^ ri ^ 6;; ^ Yt ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ rls ri * >{c ^4 rt 3?: aS 5 



90 



•m--k # 



CONTINUOUS MLD FORECAST FCR STATION NOVEMBER 



FORECAST FCR OZ 11 SEP. 

CONVECTIVE MLD IS 0.0 

FORCED ML C IS 29.02 
LAYER DEPTH IS 29.02 
CGNV/DIV CORR IS -0.44 

FORECAST MLD IS 28.57 
ACTUAL MLD IS 160.00 



FORECAST FCR OZ 12 SEP. 

CONVECTIVE MLD IS 0.0 

FORCED MLD IS 18.86 
LAYER DEPTH IS 18 .86 
CONV/DIV CORR I S -0. 20 

FORECAST MLD IS 18.66 
ACTUAL MLD IS 120.00 



FORECAST FOR OZ 13 SEP. 

CCNVECT IVE MLD IS 0.0 

FORCED MLD IS 53.71 
LAYER DEPTH IS 53.71 
CONV/DIV CORR IS -0.37 

FORECAST MLD IS 53.34 
ACTUAL MLD IS 110.00 



FORECAST FCR OZ 14 SEP. 

CONVECTIVE MLD IS 53.93 
FORCED MLD IS 55.67 
LAYER DEPTH IS 55.67 
•CONY/DIV CORR I S -C. 09 
FORECAST MLD IS 55.59 
ACTUAL MLD IS 180.00 
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CONTINUOUS MLD FORECAST FOR STATION NOVEMBER 



FORECAST FCR OZ 12 JAN. 

CONVECTIVE MLD IS 550.06 
FORCED MLD IS 132. 37 
LAYER DEPTH IS 550.06 
CON V/DI V CORR IS l .90 

FORECAST MLD IS 551.96 
ACTUAL MLD IS 530.00 



FORECAST FOR OZ 13 JAN. 

CONVECTIVE MLD IS 55L.99 
FORCED MLD IS 138.08 
LAYER DEPTH IS 551.99 
CONV/DIV CORR IS C.79 

FORECAST MLD IS 552.78 
ACTUAL MLD IS 510. CC 



FORECAST FCR OZ 14 JAN. 

CONVECTIVE MLD IS 552.61 
FOR CEC ML CIS 83.18 
LAYER DEPTH IS 552.81 
CONV/DIV CORR IS 1.61 

FORECAST MLD IS 554.42 
ACTUAL MLD IS 520.00 



FORECAST FOR OZ 15 JAN. 

CONVECTIVE MLD IS 0.0 

FORCED MLD IS 42.45 
LAYER DEPTH IS 42.45 
CONV/DIV CORR IS 0.33 

FORECAST MLD IS 42.78 
ACTUAL MLD IS 540. CO 



FORECAST FCR OZ 16 JAN. 

CONVECTIVE MLD IS 45.77 
FORCED MLD IS 82.02 
LAYER DEPTH IS 82.02 
CONV/DIV CORR IS -0.23 
FORECAST MLD IS 81.79 
ACTUAL MLD IS 440.00 
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FORECAST MLD FOR STATION NO V EMBE R ( CONT ) 



FORECAST FCR 02 17 JAN. 

CONVECTIVE MLD IS 84.14 

FGRCEC ML C IS 114.14 
LAVER DEPTH IS 114.14 
CONV/DIV CORR IS C.44 

FORECAST MLD IS 114.58 
ACTUAL MLD IS 540.00 



FORECAST FOR OZ 18 JAN. 

CGNVF.CTIVE MLD IS 0.0 

FORCED MLD IS 46.32 
LAYER DEPTH IS 46.22 
CONV/DIV CORR IS -0.57 

FORECAST MLD IS 45.75 
ACTUAL MLD IS 520.00 



FORECAST FCR OZ 19 JAN. 

CONVECTIVE MLD IS 46. 17 
FORCED MLC IS 77 .63 
LAVER DEPTH IS 77.63 
CONV/DIV CORR IS -0.29 
FORECAST MLD IS 77.35 
ACTUAL MLD I S 500. CO 



FORECAST FOR OZ 20 JAN. 

CONVECTIVE MLD IS 77.62 
FORCED ML C I S 87. 75 
LAY ER DEPTH IS 87 .75 
CONV/DIV CORR IS -0.26 
FORECAST MLD IS 87.48 
ACTUAL MLD IS 590.00 
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CONTINUOUS MLD FORECAST FOR STATION PAPA 



FORECAST FOR 0 Z 13 JAN. 

CONVECTIVE MLD IS 390.15 

FORCED MLD IS 164.92 
LAYER DEPTH I S 390. 15 
COKV/DIV CORR IS 10. L4 

FORECAST MLD IS 400.29 
ACTUAL MLD IS 360. CO 



FORECAST FOR OZ 14 JAN. 

CONVECTIVE MLD IS 0.0 

FORCED ML C IS 231. 18 
LAYER DEPTH IS 231.18 
CONY/DI V CORR IS 1.57 

FORECAST MLD IS 232.75 
ACTUAL MLD IS 9999.00 



FORECAST FOR OZ 15 JAN. 

CONVECTIVE MLD IS 0.0 

FORCED MLD IS 332.13 
LAYER DEPTH I S 332. 13 
CGNV/DIV CORR IS 4.07 

FORECAST MLD IS 336.20 
ACTUAL MLD IS 9999. CO 



FORECAST FCR OZ 16 JAN. 

CONVECTIVE MLD IS 337.39 
FORCED MLD IS 211.43 
LAYER DEPTH IS • 337 .39 

CONV/DIV CORR IS -1.25 
FORECAST MLD IS 336.14 
ACTUAL MLD IS 330.00 



FORECAST FOR OZ 17 JAN. 

CONVECTIVE MLD IS 338.42 

FORCED MLD IS 204.18 
LAYER DEPTH IS 338.42 
CONV/DIV CORR IS 1.04 

FORECAST MLD IS 339.46 
ACTUAL MLD IS 9999.00 
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FORECAST MID FOR STATION PAPA(CONT) 



FORECAST FOR OZ 18 JAN. 

CONVECTIVE ML D IS 341.38 
FORCED ML D IS L92.C8 
LAYER DEPTH IS 341.38 
CGNV/DIV CORR IS 6.34 

FORECAST KLD IS 347.92 
ACTUAL ML D IS 430.00 



FORECAST FCR OZ 19 JAN. 

CONVECTIVE MLD IS 349.43 
FORCED ML C IS 139.04 
LAYER DEPTH IS 349.43 
CONV/DIV CORR IS 7.22 

FORECAST MLD IS 356.65 
ACTUAL MLD IS 390.00 



FORECAST FOR OZ 20 JAN. 

CONVECTIVE MLD IS 358.78 
FORCED MLD IS 172.31 
LAYER DEPTH IS 358.78 
CONV/DIV CORR IS 10.02 
FORECAST MLD IS 368.80 
ACTUAL MLD IS 380.00 
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